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Abstract 



We report a new measurement of the cross section for the production of 
isolated photons, with transverse energies (Ej>) above 10 GeV and pseudo- 
rapidities \rj\ < 2.5, in pp collisions at y/s = 1.8 TeV. The results are based 
on a data sample of 107.6 pb _1 recorded during 1992-1995 with the D0 de- 
tector at the Fermilab Tevatron collider. The background, predominantly 
from jets which fragment to neutral mesons, was estimated using the longi- 
tudinal shower shape of photon candidates in the calorimeter. The measured 
cross section is in good agreement with the next-to-leading order (NLO) QCD 
calculation for Ej, > 36 GeV. 
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Direct (or prompt) photons, by which we mean those produced in a hard parton-parton 
interaction, provide a probe of the hard scattering process which minimizes confusion from 
parton fragmentation or from experimental issues related to jet identification and energy 
measurement Jl| . In high energy pp collisions the dominant mode for production of photons 
with moderate transverse energy Ej, is through the strong Compton process qg —>■ qj. The 
direct photon cross section is thus sensitive to the gluon distribution in the proton. Direct- 
photon measurements allow tests of NLO and resummed QCD calculations, phenomeno- 
logical models of gluon radiation, and studies of photon isolation and the fragmentation 
process. 

Data from previous collider measurements [@-@] have indicated an excess of photons at 
low E]n (< 25 GeV) compared with predictions of NLO QCD. This excess may originate 
in additional gluon radiation beyond that included in the QCD calculation ||, or reflect 
inadequacies in the parton distributions and fragmentation contributions ||. 

In this Letter, we present a new measurement of the cross section for production of 
isolated photons with E% > 10 GeV and pseudorapidity \r}\ < 2.5 in collisions at 
y/s = 1.8 TeV, which supersedes our previous publication [|J. (Pseudorapidity is defined 
as 7] = —In tan | where 9 is the polar angle with respect to the proton beam.) The higher 
statistical precision afforded by the increased luminosity (12.9 ± 0.7 pb _1 recorded during 
1992-1993 and 94.7 ±5.1 pb _1 recorded during 1994-1995) motivated a refined estimation 
of the backgrounds. In particular, fully-simulated jet events were used in place of single 
neutral mesons to model background. 

Photon candidates were identified in the D0 detector [|7J as isolated clusters of energy 
depositions in the uranium and liquid-argon sampling calorimeter. The calorimeter covered 
\rj\ < 4 and had electromagnetic (EM) energy resolution /Je/E m 15% /-^E (GeV) © 0.3%. 
The EM section of the calorimeter was segmented longitudinally into four layers (EM1-EM4) 
of 2, 2, 7, and 10 radiation lengths respectively, and transversely into cells in pseudorapidity 
and azimuthal angle Arj x A0 = 0.1 x 0.1 (0.05 x 0.05 at shower maximum in EM3). Drift 
chambers in front of the calorimeter were used to distinguish photons from electrons, or 
from photon conversions, by ionization measurement. 

A three-level trigger was employed during data taking. The first level used scintillation 
counters near the beam pipe to detect an inelastic interaction; the second level required that 
the EM energy in calorimeter towers of size At] x A0 = 0.2 x 0.2 be above a programmable 
threshold. The third level was a software trigger in which clusters of calorimeter cells were 
required to pass minimal criteria on shower shape. 

Offline, candidate clusters were accepted within the regions \r]\ < 0.9 (central) and 
1.6 < \tj\ < 2.5 (forward) to avoid inter-calorimeter boundaries; in the central region, clusters 
were required to be more than 1.6 cm from azimuthal boundaries of modules. The event 
vertex was required to be within 50 cm of the nominal center of the detector along the beam. 
Each candidate was required to have a shape consistent with that of a single EM shower, 
to deposit more than 96% of the energy detected in the calorimeter in the EM section, and 
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to be isolated as defined by the following requirements on the transverse energy observed 
in the annular region between 7Z = a/ At? 2 + A0 2 = 0.2 and 1Z = 0.4 around the cluster: 
j^tz<oa _ j^r.<o.2 ^ 2 G e v. The combined efficiency of these selections was estimated as a 
function of Ej, using a detailed Monte Carlo simulation of the detector || and verified with 
electrons from Z — > ee events, and found to be 0.65 ± 0.01 (0.83 ± 0.01) at Ej, = 40 GeV for 
central (forward) photons. An uncertainty of 2.5% was added in quadrature to this to allow 
for a possible dependence on instantaneous luminosity. Photon candidates were rejected 
if there were tracks within a road A6 x A0 m 0.2 x 0.2 radians between the calorimeter 
cluster and the primary vertex. The mean efficiency of this requirement was measured 
to be 0.83 ± 0.01 (0.54 ± 0.03) in the central (forward) region. The inefficiency stemmed 
mainly from photon conversions and overlaps of photons with charged tracks (either from 
the underlying event or from other pp interactions) . 

Background to the direct-photon signal comes primarily from two-photon decays of 7r° 
and 77 mesons produced in jets. While the bulk of this background is rejected by the selection 
criteria (especially the isolation requirement), substantial contamination remains, predomi- 
nantly from fluctuations in jet fragmentation, which can produce neutral mesons that carry 
most of the jet energy. For a ir° meson with Ej, ^ 10 GeV, the showers from its two-photon 
decay coalesce and mimic a single photon in the calorimeter. 

The fraction of the remaining candidates that are genuine direct photons (the purity V) 
was determined using the energy E\ deposited in the first layer (EMI) of the calorimeter. 
The decays of neutral mesons primarily responsible for background produce two nearby 
photons, and the probability that at least one of them undergoes a conversion to an e + e~ 
pair either in the cryostat of the calorimeter or the first absorber plate is roughly twice 
that for a single photon. Such showers due to meson decays therefore start earlier in the 
calorimeter than showers due to single photons, and yield larger E\ depositions for any initial 
energy. A typical distribution in our discriminant, log 10 [1 + log 10 {1 + i?i(GeV)}], is shown 
in Fig. [l]. This variable emphasized differences between direct photons and background, 
and was insensitive to noise and event pileup. A small correction, based on electrons from 
W decays, was made to bring the E\ distribution for the 1992-1993 data into agreement 
with the 1994-1995 data. The distribution in the discriminant was then fitted to the sum 
of a photon signal and jet background, both of which were obtained from Monte Carlo 
simulation. Two components of the jet background were included separately: those with 
and those without charged tracks inside the inner isolation cone (TZ = 0.2 from the photon 
candidate). This was done to minimize constraints in the fit from the (relatively poorly 
determined) tracking efficiency and from the model used for jet fragmentation. 

Direct photon and QCD jet events were generated using pythia || and then passed 
through the geant detector-simulation package, and overlaid with data acquired using a 
random trigger to model noise, pileup, underlying event, and multiple pp interactions M. 
The simulated E\ was corrected for imperfect modeling of the material in the detector. We 
assumed that the Monte Carlo energy could be parametrized as E^ c = a + /3Ex, with the 
parameters a and f3 determined from data: f3 from the W — > ev sample and a from the 
photon data. The fits to extract the purity V were performed for different values of a, and 
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Iog 10 [1 +log 10 {(1+E 1 (GeV)}] 

FIG. 1. Distribution of the discriminant variable for 21 < E2, < 26GeV central photon candi- 
dates (points with error bars), and the fitted distribution (solid curve) composed of Monte Carlo 
photons (curve labelled (a)) and jets with and without charged particles (curves labelled (c) and 
(b), respectively). The Monte Carlo curves shown here were smoothed for clarity (this was not 
done in the fitting itself). Results of these fits provide the purity V of the signal: for this bin, 
V = 0.58 ± 0.07. 



the total x 2 was minimized for all E^. 

To reduce computation time, the jet background events were preselected just after their 
generation to have highly electromagnetic jets. The background subtraction technique used 
in this analysis employs fully-simulated jet events, whereas the previous analysis modeled 
the background with isolated neutral mesons. With our increased statistics, it was found 
that individual isolated mesons could not adequately model the background. Indeed, our 
simulation shows that less than half of the background can be attributed to the presence 
of single neutral mesons within the inner isolation cones (of 1Z = 0.2). The new approach 
provided a much better description of the shower shape and isolation energy, and resulted 
in an increased estimate of the signal purity. 

Fitting was done separately for samples at central and forward regions, for each Ej, 
bin, using the package hmcmll [T(J, with the constraint that the fractions of signal and 



background were between 0.0 and 1.0. The resulting purity V and its uncertainty is shown 
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FIG. 2. The fraction of photon candidates that are direct photons as a function of EjI, for 
central and forward photons. 



in Fig. as a function of Ej,. As well as the fitting error, a systematic error was assigned to 
the use of pythia to model jets. This uncertainty was estimated by varying the multiplicity 
of neutral mesons in the core of the jet by ±10% [p| . 

The differential cross section d 2 a/dE^,dr], determined after correction for purity and 
efficiency (but not corrected for energy resolution) is shown as a function of Ej, in Fig. |3| 
and in Table |. The purity corrections were applied point by point, using the same binning 
for the cross section as for the determination of purity. The correlated errors consist of the 
quadrature sum of the uncertainties on luminosity, vertex requirements, and energy scale 
in the Monte Carlo (which are energy independent) and the model for fragmentation (large 
uncertainty at low Ej, because of the low purity in this region). The uncorrelated errors 
include the statistical uncertainty, the fitting error, and the statistical uncertainties on the 
determination of acceptance, trigger efficiency, and the efficiency of the selection criteria. 

These new measurements are ~ 20 — 30% higher than our previously published results. 
The change is well understood, and is due to the improvements in the Monte Carlo model 
used to estimate the purity, and in calculations of the acceptance and luminosity [|T^]. 

We compare the measured cross section with NLO QCD calculations using the program 
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FIG. 3. The cross section d 2 a/dE^drj for isolated photons as a function of transverse energy 
-EZ, for central and forward regions. The curves show the NLO QCD calculated cross sections. 



of Baer, Ohnemus, and Owens |T3[. This calculation includes 7 + jet, 7 + two jets, and 
two jets with bremsstrahlung in the final state. In the latter case, a jet collinear with the 
photon was created with the remaining fraction of the energy of the relevant final-state 
parton, so that the isolation cut could be modeled. For all sources of signal, the final-state 
parton energies were smeared using the measured EM and jet resolutions. The isolation 
criterion was imposed by rejecting events with a jet of Et > 2 GeV within 1Z < 0.4 of the 
photon. (Smearing photon and jet energies changed the QCD prediction by less than 4%.) 



CTEQ4M parton distributions |I4| were used in the NLO calculations, with renormalization 
and factorization scales fiR = Hf = -E™ ax , where i?™ ax is the larger of the transverse energies 
of the photon or the leading jet. If, instead, the scales fifi — fip — 2£/™ ax or _£/™ ax /2 were 
employed, the predicted cross sections changed by < 6%. 

Figure |] shows the difference between experimental and theoretical differential cross 
sections (d 2 a/dE^ dt]), divided by the theoretical values. In both central and forward regions, 
the NLO QCD predictions agree with the data for transverse energies Ej, ^ 36 GeV. At 
lower transverse energies, particularly for |?7| < 0.9, our measured cross section exceeds the 
expectation from NLO QCD, a trend consistent with previous observations at collider @-|| 
and fixed target jn| energies. Using contributions from both correlated and uncorrelated 
errors, the y 2 value for the data compared with NLO QCD is 8.9 in the central region and 
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FIG. 4. Difference between the measured differential cross section for isolated photon produc- 
tion and the prediction from NLO QCD, using CTEQ4M parton distributions. 



1.9 in the forward region, for Ej, < 36 GeV in each case (the first 4 data points). 

These data complement and extend previous measurements, and provide additional input 
for extraction of parton distributions through global fits to all data. The difference between 
the data and NLO QCD for Ej, < 36 GeV suggests that a more complete theoretical 
understanding of processes that contribute to the low-£^ behavior of the photon cross 
section is needed. 

We thank J. F. Owens for his assistance with the theoretical calculations. We thank the 
Fermilab and collaborating institution staffs for contributions to this work, and acknowledge 
support from the Department of Energy and National Science Foundation (USA), Commis- 
sariat a L'Energie Atomique (France), Ministry for Science and Technology and Ministry 
for Atomic Energy (Russia), CAPES and CNPq (Brazil), Departments of Atomic Energy 
and Science and Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry 
of Education and KOSEF (Korea), CONICET and UBACyT (Argentina), A.P. Sloan Foun- 
dation, and the Humboldt Foundation. 
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TABLE I. The predicted and measured cross sections in bins of Ej,. {E^} is the average 
photon transverse energy in each bin. The columns labelled Sajj and boc show the magnitude of 
the uncorrelated and correlated uncertainties, respectively. (The statistical error is contained in 
Sau.) 
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